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In the first part of tliis paper, we sfiow tfiat tlie semiclassical Einstein-Langevin equation, intro- 
duced in the framework of a stochastic generalization of semiclassical gravity to describe the back 
reaction of matter stress-energy fluctuations, can be formally derived from a functional method 
based on the influence functional of Feynman and Vernon. In the second part, we derive a number 
of results for background solutions of semiclassical gravity consisting of stationary and conformally 
stationary spacetimes and scalar fields in thermal equilibrium states. For these cases, fluctuation- 
dissipation relations are derived. We also show that particle creation is related to the vacuum 
stress-energy fluctuations and that it is enhanced by the presence of stochastic metric fluctuations. 
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I. INTRODUCTION 



It is generally believed that there must be a regime in which the gravitational field can be treated as a classical or 
"quasiclassical" field, but its interaction with quantum matter fields cannot be neglected. The standard approach to 
describe such a regime is the semiclassical theory of gravity based on the semiclassical Einstein equation. This is a 
generalization of the Einstein equation for a classical metric when the expectation value of the stress-energy tensor 
of quantum matter fields is the source of curvature. The semiclassical theory of gravity is mathematically consistent 
and fairly well understood, at least for linear matter fields 

One expects that semiclassical gravity could be derived as an approximation of a fundamental quantum theory of 
gravity. However, in the absence of such a fundamental theory, the scope and limits of the semiclassical theory are 
less well understood [QJ^. It has been pointed out, nevertheless, that this semiclassical theory may not be valid when 
the matter fields have important quantum stress-energy fluctuations When this is the case, the stress-energy 

^\ fluctuations may have relevant back-reaction effects on the spacetime geometry in the form of induced gravitational 
fluctuations [Q. A number of examples have been studied, both in cosmological and in flat spacetimes, where, for 
qh, some states of the matter fields, the stress-energy tensor have significant fluctuations It is thus necessary to extend 
^ • the semiclassical theory of gravity to determine the effect of such fluctuations. 

To address this problem, different approaches have been adopted. The aim of the first part of the present paper is 
• • ■ to unify two of these approaches. 

One of these approaches relies on the idea, first proposed by Hu |^ in the context of semiclassical cosmology, of 
viewing the metric field as the "system" of interest and the matter fields as being part of its "environment." This 
$^ ' approach leads naturally to the infiuence functional formalism of Feynman and Vernon [|lO|. In this formalism, the 
5^ integration of the environment variables in a path integral yields the influence functional, from which one can define 
an effective action for the dynamics of the system [ pd| p^ . This approach has been extensively used in the literature, 
not only in the framework of semiclassical cosmology []l2p ^ , p^ 25 1 , but also in the context of analogous semiclassical 
regimes in quantum mechanics |l5| , ^j26| and in quantum field theory [p^ , p7| -^0| . It is based on the observation that 



the semiclassical equation can be directly derived from the effective action of Feynman and Vernon ]12| , [18| , |19y22|j23| , |2£ 
When computing this effective action perturbatively up to quadratic order in its variables, one usually finds some 
imaginary terms which do not contribute to the semiclassical equation. The key point is then to formally identify the 
contribution of such terms in the influence functional with the characteristic functional of a Gaussian stochastic source. 
Assuming that, in the semiclassical regime, this stochastic source interacts with the system variables, equations of the 
Langevin type can be derived for these variables. However, since this approach relies on a purely formal identification, 
doubts may be raised on the physical significance of the derived equations. 

An alternative approach has been introduced in a recent paper ||3l|| . In that work, we proposed a stochastic 
semiclassical theory of gravity as a perturbative generalization of semiclassical gravity to describe the back reaction 
of the lowest order stress-energy fluctuations. The idea is in fact quite simple. One starts realizing that, for a given 
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solution of semiclassical gravity, the lowest order matter stress-energy fluctuations can be associated to a classical 
stochastic tensor field. Then, we seek an equation which incorporates in a consistent way this stochastic tensor as the 
source of linear perturbations to the semiclassical metric. The resulting equation is the semiclassical Einstein-Langevin 
equation. 

We should emphasize that, even if the metric fluctuations in this theory are classical (stochastic fluctuations), their 
origin is presumably quantum. This is so not only because these metric fluctuations are induced by the fluctuations of a 
quantum operator, but also because they are supposed to describe some remnants of the quantum gravity fluctuations 
after some mechanism for decoherence and classicalization of the metric held |32-3^. From the formal assumption 
that such a mechanism is the Gell-Mann and Hartle mechanism of environment-induced decoherence of suitably 
coarse-grained system variables [^,^, one may, in fact, derive the stochastic semiclassical theory |3^. Nevertheless, 
that derivation is of course formal, given that, due to the lack of the full quantum theory of gravity, the classicalization 
mechanism for the gravitational flcld is not understood. One expects that the stochastic semiclassical theory is valid 
when the characteristic time and space scales of variation of the metric field are well above its characteristic decoherence 
scales. In this regime, the theory can be applied to compute correlation functions of gravitational perturbations for 
points separated by scales larger than these decoherence scales. Hence, this theory may have a number of interesting 
applications in black hole physics and in cosmology, particularly in view of the problem of structure formation. Some 

examples of simple applications have already been given in Refs. [mm^. 

The purpose of the second part of the paper is to derive some general results concerning stochastic semiclassical 
gravity for stationary and conformally stationary background solutions of semiclassical gravity (for conformal matter 
fields in the latter case). We analyze two issues: the existence of a fluctuation-dissipation relation and the creation of 
particles by stochastic metric perturbations. 

Under very general conditions, a fluctuation-dissipation relation is known to exist in models of quantum mechan- 
ics, and also in some models of quantum many-body systems or quantum fields in the presence of classical fields 
p8|- p| , p2 45 1 . This is a relation between quantum fluctuations of a system in a state of thermal equilibrium and 
the dissipative properties of this system caused by classical linear perturbations on it. The idea of a fluctuation- 
dissipation relation in the theory of quantum fields in curved spacetimes and in the semiclassical back-reaction 
problem was already present in some early papers 1 46 . A fluctuation-dissipation relation has been found in some 
of the previous derivations of semiclassical Langevin-type equations ||l^ , |l9| , ^ , p4| . Some authors believe that such a 
relation should always be present and embody the physics of the back reaction of matter fields on the gravitational 
field |l|,|2|,||,||,|9|. It is also believed that noise and dissipation must be related to the creation of particles by 
stochastic metric perturbations ||,|l^-|l|,|l|,|2|j4|@ . 

In stationary and conformally stationary spacetimes (for conformal fields in the latter case), one can define a state 
of thermal equilibrium for the matter fields. When the background solution of semiclassical gravity is of one of these 
types, we can identify a dissipation kernel in the corresponding semiclassical Einstein-Langevin equation which is 
related to the fluctuations of the stochastic source by a fluctuation-dissipation relation. We also study the production 
of particles by stochastic metric perturbations to such backgrounds: we relate particle creation to the vacuum stress- 
energy fiuctuations and we show that the mean value of created particles is enhanced by the presence of metric 
fluctuations. 

The plan of the paper is the following. In Sec. |l|, we construct the stochastic semiclassical theory of gravity to 
describe the back reaction of the stress-energy fluctuations on the spacetime. In Sec. Ill, we show that the semiclassical 
Einstein-Langevin equation obtained in Sec. |l| can actually be formally derived with the functional approach. This 
connection clarifies the physical meaning of the Langevin-type equations previously derived by functional methods 
^- |l4| , p^ |2^] , since it shows that the formally introduced stochastic source is directly related to the matter stress- 
energy fiuctuations. We then use the functional approach to write the Einstein-Langevin equation in an explicit form, 
which is more suitable for specific calculations. In Sec. IV, we derive the fluctuation-dissipation relation for stationary 
and conformally stationary backgrounds and the results for particle creation by stochastic metric perturbations. 
Finally, in Sec. 0, we summarize our main conclusions. 

Throughout this paper we use the (4- 4- -I-) sign conventions and the abstract index notation of Ref. |Q, and we 
work with units in which c = H = 1. 



II. STOCHASTIC SEMICLASSICAL GRAVITY 



In this section, we construct the stochastic semiclassical theory of gravity as a perturbative extension of semiclassical 
gravity to describe the back reaction of quantum stress-energy fluctuations on the gravitational fleld. Let us begin with 
a brief overview of the semiclassical theory of gravity interacting with linear matter flelds. Let {Ai,gab) be a globally 
hyperbolic four-dimensional spacetime and consider a linear quantum fleld on it. For the sake of definiteness, we 
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will take $ as a real scalar field, but all the analysis of this section is valid for any kind of linear quantum field or 
for a set of linear independent quantum fields. Throughout this section we shall work in the Heisenberg picture. The 
field operator in this picture, is an operator-valued distribution solution of the Klein-Gordon equation, 

(□-m2-^i?)$ = 0, (2.1) 

where m is the mass, □ = VaV": with \/a being the covariant derivative associated to the metric gab, and £, is 
a dimensionless parameter coupling the field to the scalar curvature R. To indicate that the field operator is a 
fimctional of the metric gab, we will write 

The classical stress-energy tensor is obtained by functional derivation of the classical action for the field in a 
background spacetime (M, gab) with respect to the metric. This tensor is a functional Tab[g, of the metric gab and 
of the classical field $. For a real scalar field, it is 

Tab[g,^] = Va$ Vfj^"- ^gabV'^^Vc^ - ^ gab + S. (gab^ - VaVb + Gab) (2.2) 

where Gab is the Einstein tensor. The next step is to define a stress-energy tensor operator Tab[g]{x). In a naive 
way, one would replace the classical field <i> in the functional Tab[g, $] by its corresponding quantum operator <&[g]. 
However, since the field operator is well-defined only as a distribution on spacetime and this procedure involves taking 
the product of two distributions at the same spacetime point, the formal expression for Tah[g] is ill-defined and we 
need a regularization procedure. We may formally think of a regularized "operator" Tab[g]{x; fl), depending on some 
regulator $1, defined by giving a precise prescription for computing its matrix elements for physically acceptable states 
of the field. These states are assumed to be Hadamard states on the Fock space of a Hadamard vacuum state [||. 
The states may have to be regularized also in some way and the procedure may involve some analytic continuation 
in the values of the regulator. Of course, if we remove the regularization in the results for these matrix elements, we 
would obtain infinite quantities. 

Once the regularization prescription has been introduced, a renormalized and regularized stress-energy "operator" 
T^[g]{x; Q) may be defined as 

f^dgK^; n) = fab[g]{x; n) + Fg[g]{x; n) /, (2.3) 

where / is the identity operator and are some symmetric tensor counterterms, which can be written in terms 

of the regulator fl and local functionals of the metric 5c(i(2;)j^ These counterterms can and must be chosen in such 
a way that, for any pair of physically acceptable states \^) and |</?), the matrix element of the renormalized operator 
f defined by 



Wt"b\ip) ^ ^Hm jV^lf «|(p)(J7), (2.4) 

where 17„ means the "physical value" of the regulator, is finite (well defined as a distribution) and satisfies Wald's 
axioms [pl0|. Using the point-splitting or the dimensional regularization methods, these counterterms can be extracted 
from the singular part of a Schwinger-DeWitt series |^,^,g2|. The choice of these counterterms is not unique, each 
different choice is called a "renormalization scheme," and this leads to some ambiguity in the definition of the 
renormalized stress-energy tensor operator. But this ambiguity can be absorbed into the renormalized coupling 
constants appearing in the equations of motion for the gravitational field. Thus, the ambiguity is only a mathematical 
artifact of the separation of the action into a gravitational part and a matter part, but the physically relevant equations 
are in fact unique 

The semiclassical Einstein equation for the metric gab can then be written as 

{Gab[g] + Agab) ~ 2 {aAab + f3Bab)[g] = {f^b)[9], (2.5) 



^In the point-splitting regularization method, for instance, one introduces a point y in a normal neighborhood of the point x, 
so some non-local dependence on the metric is explicitly introduced in the regularized stress-energy operator and then also in 
the counterterms. Using this regularization technique, the regulator can be taken as the vector cj°'{x,y), which is the tangent 
vector at the point x to the geodesic joining x and y with length equal to the arc length along this geodesic. In this case, the 
counterterms can be written in terms of the vector a"'{x,y) and tensors which are local functionals of the metric gab{x) [S,^- 
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where (r^)[.g] is the expectation value of T^[g] in some physically acceptable state of the quantum field on the 
spacetime {A4,gab)- The notation {T^)[g] is used to indicate that this expectation value is a functional of the metric 
gcd, not only because the stress-energy tensor operator depends on the metric, but also because the state of the matter 
field depends on the spacetime (in general, such state depends on the global structure of the spacetime manifold). 
In Eq. ( ^.5| ), G, A, a and /3 are renormalized coupling constants, respectively, the Newtonian gravitational constant, 
the cosmological constant and two dimensionless coupling constants. These constants may be seen as the result of 
"dressing" the bare coupling constants in a suitably regularized version of the gravitational part of the action. 
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{R-2Ab) 



(2.6) 



where Cabcd is the Weyl tensor and the subindex b in the coupling constants means "bare." These renormalized 
coupling constants are supposed to be determined experimentally (for the specific renormalization scheme that one 
has chosen and for the characteristic scales of the physics under consideration) . The tensors Aab and Bab in Eq. (|2.5| ) 
come from the functional derivatives with respect to the metric of the terms quadratic in the curvature in Sg [g] , which 
are needed to ensure the renormalizability of the theory. These tensors are explicitly given by 
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where Rabcd is the Riemann tensor and Rab is the Ricci tensor. Note that each of the terms in Eq. ( |2.5| ) has vanishing 
divergence. Notice also that we could add a classical stress-energy tensor to the right hand side of Eq. ( ^.5[ ), if we 
had a classical matter source, but, for simplicity, we shall ignore such a term. 

As long as the gravitational field is assumed to be described by a classical Lorentzian metric gab, the semiclassical 
Einstein equation seems to be the only physically plausible dynamical equation for this metric. The reason is that, 
in classical general relativity, the metric gab couples to matter through the stress-energy tensor. For a field quantized 
on the spacetime {M,gab) and for a given state of this field, the expectation value of the renormalized stress-energy 
tensor operator is the only physically observable (up to the ambiguity mentioned above) c-number stress-energy tensor 
that we can construct. 

A solution of semiclassical gravity consists of a spacetime {A4 , gab) , a quantum field operator $ [g] satisfying Eq. ( ^.l| ) , 
and a physically acceptable state for this field (which can also be a mixed state characterized by a density 

operator), such that Eq. (2.5) is satisfied when the expectation value in the state of renormalized operator 

T^[g] is put on the right hand side. 

Let us now introduce stress-energy fluctuations. Given a solution of semiclassical gravity, the stress-energy tensor 
will in general have quantum fluctuations. To lowest order, such fluctuations are described by the bi-tensor, which 
shall be called noise kernel, deflned by 



8Nabcd{x,y) 



^hm {{iab{x)Jcd{y)})[gm, 



(2.9) 



where { , } means the anticommutator and tab{x; = Tab{x; ^l) — {Tab{x)) (fl) . Note that we have defined this noise 
kernel in terms of the unrenormalized "operator" Tab[g]ix;Q). For a linear quantum field, this can be done because 
the ultraviolet singular behavior of {'Tab{x)Tcd{y)) i^) is the same as that of (Tab{x ) )(Q)( Tcd{'y ))(^), so Nabcd{x,y) 
is free of ultraviolet divergencies. One can trivially see from the substitution of (|2.3| ) into ( |2.9[) that we can replace 
Tab[g]{x;Q) by the renormalized operator T^[g]{x), and omit the limit fi— >r2p, in the last expression. The result is 
obviously independent of the renormalization scheme that one chooses to define T^. 

As a perturbative correction t o se miclassical gravity, we want now to introduce an equation in which the stress- 
energy fiuctuations described by (2^) are the source of classical gravitational fluctuations. Thus, we assume that the 
gravitational field is described by gab + hab, where hab is a linear perturbation to the background metric jiab, ^solution 
of Eq. (|2.5| ). The renormalized stress-energy operator and the state of the quantum field may be denoted by T^^[g + h] 
and 1-0) [5 + h], respectively, and {T^)[g + h] is the corresponding expectation value. 

Let us introduce a Gaussian stochastic tensor field ^ab defined by the following correlators: 
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{£.ab{x)}c = 0, {£.ab{x)^cd{y))c = Nabcd{x , y) , 



(2.10) 



where ( )c means statistical average. In general, the two-point correlation function of a stochastic tensor field S^ab 
must be a symmetric, in the sense that {£,abix)S,cd{y))c = {£,cdiy)£.ab{x))c, and positive semi-definite real bi-tensor 
field. Since the renormalized operator is self-adjoint, it is easy to see from the definition (2.9) that Nabcd{x,y) 



satisfies all these conditions. Therefore, the relations ( |2.10| ), with the cumulants of higher order taken to be zero, do 
truly characterize a stochastic tensor field S,ab- The simplest equation which can incorporate in a consistent way the 
stress-energy fluctuations described by Nabcd{x, y) as the source of classical metric fluctuations is 



SttG 



Gab[9 + h]+K {gab + hab) - 2 {aAab + PBab) [g + h] = (fil) [g + h]+ 2U, 



(2.11) 



which must be understood as a dynamical equation for hab to linear order. Eq. (2.11) is the semiclassical Einstein- 
Langevin equation, which gives a first order correction to semiclassical gravity. One could also seek equations describ- 
ing higher order corrections, which would involve higher order stress-energy fluctuations, but, for simplicity, we shall 

stick to the lowest order. 

In order to check the consistency of Eq. (2.11), note that the term ^afc does not depend on hcd, since it is completely 



determined from the solution of semiclassical gravity by the correlators ( 2.10| ). Even so, this term must be considered 
as of first order in perturbation theory around semiclassical gravity. As shown in Ref . , ^ab is covariantly conserved 
up to first order in this perturbation theory, in the sense that V^Cafe behaves deterministically as the zero vector field 
on A4 (v° is the covariant derivative assoc iated to the background metric gab)- It is thus consistent to include the 
term ^ab in the right hand side of Eq. (^.ll[ ) . 

It was also shown in Ref. that for a conformal field, i.e., a field whose classical action is conformally invariant 
{e.g., a massless conformally coupled scalar field), is "traceless" up to first order in perturbation theory, since 
9°'^£,ab behaves determi nistic ally as a vanishing scalar. Hence, in the case of a conformal matter field, the trace of the 
right hand si de of Eq. (2.11) comes only from the trace anomaly. 

Since Eq. ( 2.11 ) is a linear stochastic equation for hab with an inhoniogeneous term S,ab, a solution can be formally 
written as a functional h ab\S} . Such a solution can be characterized by the whole family of its correlation functions. 
From the average of Eq. ( 2.11 ), the average of the metric, gab + {hab)c, must be a solution of the semiclassical Einstein 
equation linearized around gab- The fluctuations of the metric around this average can be described by the moments 
of order higher than one of the stochastic field ^^^[C] = hab[£,] — {hab)c- 



Finally, for the solutions of Eq. (2.11) we have the gauge freedom hab — ^ h'^f, = ^"-b + VaC& + VbCa, where is 
any stochastic vector field on A4 which is a functional of ^cd, and Ca = gabC^- Note that the tensors which appear in 
Eq. (2.11) transform as Rab[g + h'] = Rab[g + h] + £^Rab[g] (to line ar ord er in the pert urba tions), where £^ is the Lie 
derivative with respect to C"^- If we substitute hab by h'^b in Eq. (2.11), we get Eq. (2.11) plus the Lie derivative of 
a combination of the tensors which appear in Eq. (p.5|). This last tensorial combination vanishes when Eq. (2.5) is 



satisfied. Thus, it is necessary t hat th e set {M,gab, ^[g], Wid]) be a solution of semiclassical gravity to ensure that 
the Einstein-Langevin equation (2.11) is gauge invariant. 



III. DERIVATION FROM AN INFLUENCE ACTION 



The purpose of this section is to derive the semiclassical Einstein-Langevin equation (2.11) by a method based on 



functional techniques. The same method has been in fact used in the literature to derive Langevin-type equations in 



the context of semiclassical cosmology |12|- |14| , [19 -25 and of analogous semiclassical regimes for systems of quantum 
mechanics ]r5| , p^ , p6[ and of quantum field theory ]Tq , p7| -pOt . Using these functional techniques, we also work out the 
Einstein-Langevin equation more explicitly, in a form more suitable for specific calculations. Here, we consider again 
the simplest case of a linear real scalar field $. 

These functional techniques are based on the closed time path (CTP) functional formalism, due to Schwinger and 
Keldysh This formalism is designed to obtain expectation values of field operators in a direct way and it 

is suited to derive dynamical equations for expectation values; see Refs. [HmQ for detailed reviews. In our case, 
this formalism will be useful to obtain an expression for the expectation value {T'^^)[g + h] as an expansion in the 
metric perturbation. When the full quantum system consists of a distinguished subsystem (the "system" of interest) 
interacting with an environment (the remaining degrees of freedom), the CTP functional formalism turns out to 
be related |T^ , p^ , p^ , p3| , p8| , p9p7| ] to the infiuence functional formalism of Feynman and Vernon ||l^. In this latter 
formalism, the integration of the environment variables in a CTP path integral yields the influence functional, from 
which one can define an effective action for the dynamics of the system [|TT|-p^. Applying this influence functional 



formalism to our problem, the semiclassical Einstein-Langevin equation will be formally derived in subsection [II B 
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In our case, we consider the metric field gab{x) as the "system" degrees of freedom, and the scalar field ^{x) and 
also some "high- momentum" gravitational modes as the "environment" variables. Unfortunately, since the form 
of a complete quantum theory of gravity interacting with matter is unknown, we do not know what these "high- 
momentum" gravitational modes are. Such a fundamental quantum theory might not even be a field theory, in which 
case the metric and scalar fields would not be fundamental objects [Q. Thus, in this case, we cannot attempt to 
evaluate the influence action of Feynman and Vernon starting from the fundamental quantum theory and performing 
the path integrations in the environment variables. Instead, we introduce the influence action for an effective quantum 
field theory of gravity and matter [S^Jl^ , in which such "high- momentum" gravitational modes are assumed to have 
been already "integrated out." Adopting the usual procedure of effective field theories |^,^, one has to take the 
effective action for the metric and the scalar field of the most general local form compatible with general covariance: 
S[g,^] = Sg[g] + S^[g,^] + ■■■, where Sg[g] is given by @, 

5„,[.g,$] =-i Jd\^[g'''da<S>db^+{m^ + ^R)<f^] , (3.1) 

and the dots stand for terms of order higher than two in the curvature and in the number of derivatives of the scalar 
field [because of the Gauss-Bonnet theorem in f our spacetime dimensions, no further terms of second order in the 



curvature are needed in the gravitational action (2.6)]. In this paper, we shall neglect the higher order terms as well 
as self-interaction terms for the scalar field. The second order terms are necessary to renormalize one-loop ultraviolet 
divergencies of the scalar field stress tensor. Since is a globally hyperbolic manifold, we can foliate it by a family 
of inconstant Cauchy hypersurfaces St. We denote by x the coordinates on each of these hypersurfaces, and by ti 
and t f some initial and final times, respectively. The integration domain for the action terms must be understood as 
a compact region U of the manifold A^, bounded by the hypersurfaces St- and St^,. 

Assuming the form (3.1) for the effective action which couples the scalar and the metric fields, we can now introduce 
the corresponding infiuence functional. This is a functional of two copies of the metric field that we denote by g^j, 
and g~ij. Let us assume that, in the quantum effective theory, the state of the full system (the scalar and the metric 
fields) in the Schrodinger picture at the initial time t = ti can be described by a factorizable density operator, i.e., a 
density operator which can be written as the tensor product of two operators on the Hilbert spaces of the metric and 
of the scalar field. Let p^{ti) be the density operator describing the initial state of the scalar field. If we c onsider the 
theory of a scalar field quantized in the classical background spacetime {M,gab) through the action (|3.lD, a state in 
the Heisenberg picture described by a density operator p[g] corresponds to this state. Let {|(/3(x))^} be the basis of 
eigenstates of the scalar fleld operator $^(x) in the Schrodinger picture: <I>^(x) |(^)^ = 'fii^) \ f)^- The matrix elements 
of p^{ti) in this basis will be written as Pi[ip, 0\ = ^{(p\ p^{ti) \0)^. We can now introduce the influence functional as 
the following path integral over two copies of the scalar fleld: 

^iF[g+,g-] ^ l^i^+l p^[^+{U),^^{U)] S[<P+{tf)-<i>^{tj)] e'(^^™[9^*+l-5™[9-^*-l). (3.2) 

The above double path integral can be rewritten as a closed time path (CTP) integral, namely, as an integral over 
a single copy of field paths with two different time branches, one going forward in time from U to tf, and the other 
going backward in time from tf to ti. From this influence functional, the influence action, SiY[g^ , g^], and the 
effective action of Feynman and Vernon, Scs[g^,g^], are deflned by JFiF[5+,g~] = e*'^^'''^^'^ ' and S'off[5^,5^] = 
Sglg"^] - Sci\g-] + SiY\g+,g^]. 

Expression (^.2[) is ill-deflned, it must be regularized to get a meaningful infl uen ce functional. We shall assume that 



we can use dimensional regularization, that is, that we can give sense to Eq. (3.2) by dimensional continuation of all 
the quantities that appear in this expression. We should point out, nevertheless, that for this regularization method 
to work one must be able to perform an analytic continuation to Riemmanian signature |]59| . Thus, we substitute the 



action Sm in ( |3.2| ) by some generalization to n spacetime dimensions, which may be chosen as 

Sr,, [g, $„] = - i jd^'x ^ [g'^'^da^ndb^n + (m^ + CR) ^l] , (3.3) 

where we use a notation in which a subindex n is attached to these quantities that have different physical dimensions 
than the corresponding physical quantities in four dimensions. A quantity with the subindex n can always be associated 
to another without this subindex by mea ns o f a mass scale p; thus, for the scalar fleld <&„ ^/j,*-""^-*/^ $. 



We also need to substitute the action (2.6) by some suitable generalization to n spacetime dimensions. We take 



Sg[g]=p''-' d^'x 



1 2 

(i? — 2Ab) + — as {RabcdR'^ ^ — RabR'^ ) + PbR^ 
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(3.4) 



6 



By the Gau ss-Bonnet theorem, this action gives for n = 4 the same equations of motion as the action ( ^.6| ) . The 
form of (3^) is suggested by the Schwinger-DeWitt a naly sis of t he u ltraviolet divergencies in the matter stress-energy 
tensor using dimensional regularization js^ . Using (3.3) and (3.4), one can w rite the ef fecti ve action of Feynman 
and Vernon, Scff[g^ , g~], in dimensional regularization. Since the action terms (3.3) and (3.4) contain second order 
derivatives of the metric, one should also add some boundary terms ||5^,|l^. The effect of these terms is to cancel out 
the boundary terms which appear when taking variations of Scs[g^ , g~] keeping the value of g'^j^ and g^i^ fixed on the 
boundary oiU. Alternatively, in or der t o obtain the equations of motion for the metric in the semiclassical regime, 
we can work with the action terms (3.3) and (3.4) (without boundary terms) and neglect all boundary terms when 
taking variations with respect to g^i^. From now on, all the functional derivatives with respect to the metric will be 
understood in this sense. 



A. The semiclassical Einstein equation in dimensional regularization 



From the action (3^), we can define the stress-energy tensor functional in the usual way 

2 dSm[g,^n] 



\/-g{x) Sgab{x) 



(3.5) 



which yields (2.2). Working in the Heisenberg picture, we can now formally introduce the regularized stress-energy 
tensor operator as 

f^'igj^T^'ig^^M f^^Lgl^M-^^-'^rriff], (3.6) 

where (a;) is the field operator, which satisfies the Klein-Gordon equation ( |2.l[ ) in n spacetime dimensions, and 
where we use a symmetrical ordering (Weyl ordering) prescription for the operators. Using the Klein-Gordon equation, 
the stress-energy operator can be written as 



where V^^lg] is the differential operator 



V 



ah 



From the definitions (3.2), (3.5) and (3.6), one can see that 

2 5Sw[g^,g~] 



{T:\x))[g] = 



V~g(x) 5g^^{x) 



(3.7) 



(3.8) 



(3.9) 



a^=g =a 



where the expectation value is taken in the n-dimensional spacetime generalization of the state described by p{g\. 
Therefore, differentiating S'off [g"*", 5~] = Sg[g^\ — Sg[g~] + S'lp [(?''', 5~] with respect to g^^, and then setting g'^^^ = 
9ab ~ 9a-b, we get the semiclassical Einstein equation in dimensional regularization: 



^— [C^^g] + Asg"") - {\ asD'^'' + 2PbbA [g] = ^^-(^-'HT:')[g] 



SttGb 



(3.10) 



where 



D"' = fd'^x {RcdefR''''^ - RcdR''') = \ g''\RcdesR''''f - RcdR'"' + OR) - 2R^''''R'' 



cde 



2R'"''"^R^d + AR'^'rJ' -iUR"^ + v" V''^, (3.11) 



and S"'' is defined as in (2.8) but for n spacetime dimensions, although its explicit expression in terms of the metric 
and curvature tensors is the same. When n — A^ one has that D"''' — (3/2)^"'', where A'^^ is the tensor defined in (p.7|). 
From equation ( ^.10 ), renormalizing the coupling constants to eliminate t he " divergencies" in fi~^'^~'^\T^'')[g], and 
taking the limit n^4, we get the physical semiclassical Einstein equation 
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B. A formal derivation of the semiclassical Einstein-Langevin equation 



In the spirit of the previ ous se ction, we now seek a dynamical equation for a linear perturbation hat to a semiclassical 
metric gab, solution of Eq. (3.10) in rt spacetime dimensions. From the result of the previous subsection, if such equation 
were simply a linearized semiclassical Einstein equat ion , it could be obtained from an expansion of the effective action 
Ses[g + ,g + h~]. In particular, since, from Eq. (3.9), we have that 



{f^\x))[g + h] = 



5Sw [g + h+, g + h' 



^-Act{g + h){x) 6h+^{x) 



(3.12) 



the expansion of {ff^'')[g+h] to linear order in hat can be obtained from an expansion of the influence action Sip[g + 
, g + h^] up to second order in h^j^. 

To perform the expansion of the influence action, we have to compute the first and second order functional derivatives 
of SiF[g~^ , g~] and then set gab~ 9ab~ d^^b- If we do so using the path integral representation (3.2), we can interpret 
these derivatives as expectation values of operators. The relevant second order derivatives are 



1 



5^SiF[g+,g- 



1 5^SiA9+,g-] 



= -Hi: 



^/-g{x)^-g{y) 5g+^[x)5g^^[y) 

where 

Nt^%]{x,y)^\{{if{x),i'i,\y)])[g], 

'[f^\x),f^%)])[g], 



"'"[gMx.y) - iff ='^[5](x,y) + z7V,f=''[5](:r,y), 
<^'^[g](x,y)-iiVr^[5](x,y), 



(3.13) 



H:^'%]{x,y) 
K'-'mx^y) 



1 



Im 



(T\f^\x)f:\y)))[gl 

5'^S„i[g,^r. 



\/-9{x)\/-g{y) \^9ab{x)5gcd{y) 




with £f = T^^ — {T^''), and using again a Weyl ordering prescription for the operators in the last of these expressions. 
Here, [ , ] means the commutator, and we use the symbol T* to denote that, first, we have to time order the field 
operators $„ and the n ap ply the derivative operators which appear in each term of the product T°'''(x)T'^'^{y), where 
T"'' is the functional (2^). For instance. 



T* vf^n (x) (x) V/ {y)vf^n (y) 



= lim v"\2'\2X^T(4„(a:i)<l„(a;2)l'„(x3)<I'„(a;4) 



(3.15) 



where T is the usual time ordering. This T* "time ordering" arises because we have path integrals containing 
products of derivatives of the field, which can be ex pressed as derivatives of the path integrals which do no t contain 
such derivatives. Notice, from the definitions (3.14), that all the kernels which appear in expressions (3.13) are real 
and that H^'^ is al so fre e of ultraviolet divergencies in the limit n^A. 

From (3.13) and (3.14), it is clear that the imaginary part of the influence action, which does not contribute to 
the semiclassical Einstein equation ( ^.10 ) beca use th e expectation value of T^^[g] is real, contains information on the 
fluctuations of this operator. From (|3^ ) and ( |3.13| ), taking into account that SYp[g,g] = and that SYp[g~ ,g~^] — 
— S*-p[g^ t9~]i we can w rite t he expansion for the influence action Swlg-Vh^ ,g+h^] around a background metric gab 
in terms of the kernels (|3.14 ). Taking into account that these kernels satisfy the symmetry relations 



H,^^^\x, y) = H,^^'^\y, x), H^^^^^ix, y) = -Hf^^\y, x), iff =^(x, y) = K^J^\y, x), 
and introducing a new kernel 

H^'^'^ix, y) ^ H^^^'ix, y) + Hf^^\x, y), 
this expansion can be flnally written as 



(3.16) 



(3.17) 
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2 

i 

+ 2 

where we have used the notation 



I f d-xd-y^^)^/^) [K,{x)] {Ht'%]{x,y) + K:!>^%]{x,y)) {KAv)} 



{hab} = h+i^ + h^^ 



(3.18) 



(3.19) 



We are now in the position to carry out the formal derivation of the semiclassical Einstein-Langevin equation. The 
procedure is well known p2|- p^ , p^ - p5|p^p7| , p^ , p^ , p7[ -p0|| , it consists of deriving a new "improved" effective action 
using the the following identity: 



-i J d^xdTy ^ -g{x)^ -g(y) [hab(x)] N^'""^{x,y) [/led (a)] _ / ] K ] 6* / ''"^ \/-9(x) i,'^'' (x) [h^bix)] 



(3.20) 



where 'P[£,n] is the probability distribution functional of a Gaussian stochastic tensor characterized by the corre- 
lators 



(Cf(^)>c=0, 



(3.21) 



with N^J"^"^ given in (3.14), and where the path integration measure is assumed to be a s calar under diffeomorphisms 
of {Ai,gab)- The above identity follows from the identification of the right hand side of ( 3.20 ) with the characteristic 
functional for the stochastic field In fact, by differentiation of this expression with respect to \h ab], i t can be 
checked that this is the characteristic functional of a stochastic field characterized by the correlators ( 3.21 ). When 
N^'"^'^{x,y) is strictly positive definite, the probability distribution functional for is explicitly given by 



where N~l^^j^[g]{x,y) is the inverse of N^'^'^'^[g]{x,y) defined by 

1 /Xaxb , xaxb\ ^""i^-y) 



(3.22) 



d-z V^)Nf^f{x, z)N-'{z, y)^- {5'^J\ + 



(3.23) 



Using the identity (3.20), we can write the modulus of the influence functional in the approximation (3.18) as 



I ^IF [g + h'^ ,9 + h ] I 



ImSiF[s+/i+.3+'j"] ^ / i J d"x^-g(x)i^''(x) [hatix)] 



(3.24) 



where ( )c mea ns sta tistical average over the stochastic tensor Thus, the effect of the imaginary part of the 
influence action ( 3.18|) on the corresponding influence functional is equivalent to the averaged effect of the s tochastic 
source coupled linearly to the perturbations ft.^^. The influence functional, in the approximation ( |3.18 ), can be 
written as a statistical average over 



TlF[g + h+,g + h-] = /e^-^fFlh+.h-:g;i. 



with 



A1§[h+,h-;g;^n]=ReSiF[g+h+,g + h-]+ d''x ^-g(x) ^ {x) [hab{x)] + Oih") 



(3.25) 



(3.26) 



where ReS'ip can be read from the expansion (3.18). Note that the stochastic term in this action contains the 
information of the imaginary part of 5if. Introducing a new "improved" effective action 



Acslh^ ,h \g\£,n] = Sg[g 



S,[g + h-]+A'i§[h+,h-;g;^„l 



(3.27) 
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where Sg[g + h^] has to be expanded up to second order in the perturbations ft.^^, the equation of motion for hab can 
be derived as 



SAcs[h+,h ;g;^ri 



^-detig + h){x) Sh+f,{x) 



(3.28) 



h+=h'- 



From ( 3.12| ), taking into account that only the real part of the influence action contributes to the expectation value 
of the stress-energy tensor, we get, to linear order in hab. 



SttG, 



■aBD''*' + 2l3BB''*'][g+h] = ^l 



= ,,-("-4) 



(f,f)[5 + ft] + 2/i-("-4)ef, (3.29) 



where h = g'^'^g hcd, that is, g — h + 0{h ) is the inverse of the metric gab + hab- This last equation is the 
semiclassical Einstein-Langevin equation in dimensional regularization. As we h ave pointed out in section^, the two- 
point correlation function of the stochastic source in this equation [see Eq. (3.21)], given by the noise kernel defin ed in 
( |3.14 ), is free of ultraviolet divergencies in the limit n^4. Therefore, in the Einstein-Langevin equation ( |3.29| ), one 
can perform exactly the same renormalization procedure as for the semiclassical Einstein equation ( |3.1C ). After this, 
Eq. ( 3.2£ ) will yield the physical semiclassical Einstein-Langevin equation ( 2.11 ). The derivation pres ented in this 
paper clarifies the p hysical meaning of the stochastic source formally introduced in the effective action (3.26) b y the 
ident ificat ion ( 3.24 ), since it links its two-point correlation function to the stress-energy fluctuations by Eqs. ( 3.21 ) 
and ( |3.14| ). 

There is also a connection between the equations obtained by this formal functional method and the equations 
derived from the (in general, also formal) assumption that decoherence and classicalization of suitably coarse-grained 
system variables is achieved through the mechanism proposed by Gell-Mann and Hartle |Q in the consistent histories 
formulation of a quantum theory. This last approach allows to evaluate the probability distribution associated to such 
decoherent variables, given by the diagonal elements of a decoherence functional, and, under some approximations, 
to derive effective quasiclassical equations of motion for them. These effective equations of motion can be shown to 
coincide |Q with the semiclassical equations for the background and the Langevin-type equations for perturbations 
obtained from the above functional method. Taking this connection into account, we can also conclude that, if one 
formally assumes that the Gell-Mann and Hartle mechanism works for the metric field, one is lead to the semiclassical 
Einstein equation and the semiclassical Einstein-Langevin equation for the background metric and for the metric 
perturbations, respectively p7| ]. 

We end this subsection with some comments on the relation between the semiclassical Einstein-Langevin equa- 
tion ( 2.11) and the Langevin-type equations for stochastic metric perturbations recent ly d erived in the literature 
U lllgpll]. In these previous derivations, one starts with the influence functional ( |3.2[ ), with the state of the 
scalar field assumed to be an "in" vacuum or an "in" thermal state, and computes explicitly the expansion for the 
corresponding influence action around a specific metric background. One then applies the above formal method to 
derive a Langevin equation for the perturbations to this background. However, most of these derivations start with 
a "mini-superspace" model and, thus, the metric perturbations are assumed from the beginning to have a restrictive 
form. In those cases, the derived Langevin equations do not correspond exactly to our equation, Eq. (2.11), but 
to a "reduced" version of this equation, in which only some components of the noise kernel in Eq. (2.10) (or some 
particular combinations of them) influence the dynamics of the metric perturbations. Only those equations which 
have been derived starting from a completely general form for the m etric perturbations [ p^|p^j2^ , |2^ are actually 
particular cases of the semiclassical Einstein-Langevin equation (2.11). Note, however, that the stochastic equation 



derived in Refs. |23p4|] do not correspond exactly to Eq. (2.11), since the background (Minkowski spacetime and a 



scalar field in a thermal state) is not a solution of semiclassical gravity. In this case, for the reasons explained in 
Sec. 0, the equation for the metric perturbations is not gauge invariant. 



C. Explicit linear form of the Einstein-Langevin equation 



We can write Eq. ( |3.29| ) in a more explicit form by working out the expansion of (T!^'') [g + Zi] up to linear order in 
the perturbation hab- From Eq. ( ^.12| ), we see that this expansion can be easily obtained from (3. IS). Noting, from 
(|3.14D, that 



1 



Kf^%]{x,y)^--{f:\x))[g] 



ST-'^[g,<S>^](x) 



2 V-9iy) \ ^9cd{y) 



(3.30) 
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we get 



{f,fix))[g + h] ^ (f,f (x))[.9] + {n^^^''[g;h]{x))[g] - 2 d^y ,/^)H^''^%]ix,y)hM + Oih^) 



(3.31) 

where the operator T^^°'^ is defined from the term of first order in the expansion of T"'^[g + h, $„] as 

T"^[g + /i, a>„] = T^%, $„] + T<^'"^[5, $„; h] + 0(/i2), f '^"'^[g; /i] ee T'^>'^''[5, I>„[5]; /i], (3.32) 

using, as always, a Weyl o rderi ng prescription for the operators in the last definition. Note that the third term on 
the right hand side of Eq. (3.31) is a consequence of the dependence on hcd of the field operator ^n[g + h] and of the 

density operato r p\g + h] . 

Substituting ( [3.31 ) into (3.29), and taking i nto a ccount that gab satisfies the semiclassical Einstein equation (3.10), 
we can write the Einstein-Langevin equation (3.29) as 

(G<^"^''[g; h]{x) ~ Kb h'^\x)) ~ \ aB/?<^"^'[ff; h]{x) - 2PbB''^^%- h]{x) 



(3.33) 

In the last equation we have used the superindex (i) to denote the terms of first order in the expansion in hab of the 
tensors G°''[g + h], D''''[g + h] and B''''[g + h]. Thus, for instance, G"''[g + h] = G''''[g] + G'^'^^ig; h] + Oih"^). The 
explicit expressions for the tensors G^'-' jg; h], D'^-^°'^[g; h] and h] can be found in the Appendix of Ref. ^7\ , 

and r<^'°^[(7, <!>„; h] is given in Appendix W. From r<^'°^[(7, /i], we can write an explicit expression for the operator 
T//'"''. In fact, using the Klein-Gordon equation, and expressions ( |3.7D and (3.8) for the stress-energy operator, we 
have 



(3.34) 



where T°'^[g\ h] is the differential operator 



(3.35) 



It is understood that indices are raised with the background inverse metric g""^ and that all the covariant derivatives 
are associated to the metric gab- Substituting expression ( ^.34 ) into Eq. (3.33), and using the semiclassical equation 
( |3.10 ) to get an expression for /i^*^"^*^(T',"'')[(7], we can finally write the semiclassical Einstein-Langevin equation in 
dimensional regularization as 



%-kGb 
4 



G 



(l)afc_ 



g''''G'"^Kd + G'^'h"^ + G'^h"^ + As ( /i' 



bcj^a 



(x) 



as 



jj(l)ab _ 1 g^bjjcdf^^^ _^ j^acf^b _^ J^bcf^a^ _ 2(3b (^B^'^"'' g'''' B"^ h^d + B'^'^hl + B''' h"^^ (x) 



($^(a;))[g] +2yd"yv/^/.-("-4)<^'^[5](a;,y)/ied(y) = 2^-(--^^C\^), (3-36) 

where the tensors G°^, D"'' and are computed from the semiclassical metric gab, and where we have omitted 
the functional dependence on gab and hab in G*^'°^, D^'-^"-'', iJfi'"'' and J-""'' to simplify the notation. Notice that, 
in Eq. (3.36), all the ultraviolet divergencies in the limit n— >4, which must be removed by renormalization of the 
coupling constants, are in ($^(a;)) and the symmetric part H^^'^'^{x,y) of the kernel H^'"^'^{x,y), whereas the kernels 
N^'"^'^{x,y) and H^'"^'^{x,y) are free of ultraviolet divergencies. These two last kernels can be written in terms of 
F:'''''[g]{x,y) = {if{x)i';f{y))[g] as 



Nf^%]{x,y)) = ^ReF:''^''[g]{x,y), 



H:^^^%]{x,y) = ^lmF:''^%]{x,y), 



(3.37) 



where we have used that 2 {i^''{x) i^^{y)) = {{if!'{x), i^f{y)}) + { [i^ix), iff{y)]), and the fact that the first term on 
the right hand side of this identity is real, whereas the second one is pure imaginary. Once we perform the renormal- 
ization procedure in Eq. ( 3.36| ), setting n = 4 will yield the physical semiclassical Einstein-Langevin equation. Note 
that, due to the presence of the kernel H^^'^'^{x, y), this equation will be usually non-local in the metric perturbation. 
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D. The kernels for a vacuum state 



We conclude this sect ion by con sideri ng the case in which the expectation values that appear in the Einstein- 
Langevin equation ( 3.36| ) [see Eqs. ( 3.14| )] are taken in a vacuum state |0) (for a field quantized on {M,gab) in the 
Heisenberg picture), such as, for instance, an "in" vacuum. In this case we can go further and write these expectation 
values in terms of the Wightman and Feynman functions, defined as 

(3.38) 



G+{x,y) = (0| <i>„(x)$„(y) |0)[5], ^G,Jx,y) = (0| T[<^^ix)^^{y)j Mg]. 

These expressions for the kernels in the Einstein-Langevin equation will be very useful for explicit calculations. To 
simplify the notation, we omit the functional dependence on the semiclassical metric gab, which will be understood 
in all the expr essions below. 

From ( 3.37 ), we see that the kernels N^'"^'^{x,y) and H^^'^'^{x,y) are the real and imaginary parts, respectively, of 
F^'""^{x, y) = (0| f^^{x) f^'^iy) |0)- (0| f^''{x) |0)(0| f^'^iy) |0). Since, from (|3), we can write the operator f^'' as a 
sum of terms of the form |ylj,$„(a;), Bx^nix)^, where Ax and Bx are some differential operators, we can express 
pabcd^^^ y) in terms of the Wightman function using 

(^{Ax^n{x),Bx^n{x)}{Cy^n{y),'Dy^n{y)})-(^{Ax^n{x),Bxk^{x)}) (^\Cy^n{y) ,Vy^n{y)} 

= AAxCyG+ix, y) BxVyG+{x, y) + 4^,P,G+(x, y) BxCyG^ix, y), (3.39) 

where Cx and P^, are also some differential operators and where the expectation values are taken in the vacuum |0). 
This identity can be easily proved using Wick's theorem or by writing the operator <I>„(a;) in terms of the creation 
and annihilation operators of the Fock representation corresponding to the vacuum |0). Using a Schwinger-DeWitt 
expansion for the Wightman function G^(x, y), one can actually see that the two terms on the right hand side of the 
last expression are free of ultraviolet divergencies in the limit n—^A. Finally, we find 

F^'^'ix^y) = ^^^^G+Ax.,y) ^'^^G+ix^y) + ^^^''G+{x,y) V^'v;G+(x,y) 

+ 2V:'{^^G+ix,y)^'GtAx,y))+2V;-'{^'^G+ix,y)^''G+ix,y))+2V:'v:%Gr (3.40) 



where is the differential operator (B.8). From this expression and the relations (3.37), we get expressions for the 
kernels N^'"^'^(x,y) and H^^'^{x,y) in terms of the Wightman function G^{x,y). 

The kernel H^^'"^{x, y), defined in ( ^4| ), can be written in terms of the Feynman function noting that, from Wick's 
theorem. 



(3.41) 



lm(T*(^{Ax^n{x),Bx^n{x)}{Cy^niy),Vy^n{y)} 

AxCyGpJx, y) BxVyGpJx, y) + AxVyGp^x, y) BxCyGp^x, y) 



-4Im 



where, again. Ax, Bx, Cx and 'Dx are real differential operators and the expectation value is in the vacuum |0). The 
kernel H^^'^'^{x, y) is then obtained by adding up the contribution of all the differential operators which appear in the 
product T°-^{x)T'^'^{y), where T"'' is the functional (2.2). After a long calculation, we get 

HsT''i=^,y) = -\lm[vfv;G,Jx,y) y^\fG,Sx,y) + y/X^'^Jx, y) V.^^G^Jx,?;) 

~9"\x) \/f\/AG^Ax,y) \j,-\jfG^Ax,y)-g"\y) X/f \// G^Jx , y) \/J' \/yG^Jx,y) 

+ ^5"''(a:)g^'(y) v:v/G,„{x,y) VeV/G,Jx,y)+IC:\2^^G,Jx,y) v/G,„{x,y) 

-g'\v) V;G,Ax,y) SJ,yG,,,{x,y)) +lCf{2sj:G,Ax,y) vJ'G,Ax,y) 

-g-\x) v/a„(x,y) V/a„(x,y)) + 21CflCf {Gl{x,y))\ , (3.42) 



where /Cf is the differential operator 



^ i {g''\x)U^ - v/\2' + G'^\x)) - - m'g'^\x). 



(3.43) 
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An alternative expression for {x,y), which is more similar to expression (3.40), can be obtained taking into 



account that GF^{x,y) is a Green function of the Klein-Gordon equation in n spacetime dimensions, which satisfies 



(a -m'- m^)) a„(x, y) = L-^ , (3.44) 



and using that in dimensional regularization [5"(x— y)]^ = 0. Finally, note that, in the vacuum |0), the term (<i>^(x)) 



in equation ( 3.36 ) can also be written as ($^(a;)) = iGp^{x,x) = G^{x,x). 

It is worth noting that, when the points x and y are spacelike separated, ^n{x) and $«(?/) commute and, thus, 
G^{x,y) —iGp^{x,y) = (1/2)(0| {$„(a;), $„(y)} |0), which is real. Hence, from the above expressions, we have that 
H^^'^{x,y) = H^^^^lx^y) = 0. This fact is not surprising since, from the causality of the expectation value of the 
stress-energy operator, we know that the non-local dependence on the metric perturbation in the Einstein-Langevin 
equation must be causal. 

IV. FLUCTUATIONS IN STATIONARY AND CONFORMALLY STATIONARY BACKGROUNDS 

In this section, we derive a number of results concerning the stochastic semiclassical theory of gravity for two 
classes of background solutions of semiclassical gravity. The first class consists of a stationary spacetime and a scalar 
field in thermal equilibrium or in its vacuum state. In the second class, the spacetime is conformally stationary, the 
scalar field is massless and conformally coup led, and its state is the conformal vacuum or a thermal state built on 
the conformal vacuum. In subsections IV A and IV B, we identify a kernel in the corresponding Eins tein-Langevin 
equations which is related to the noise kernel by a fluctuation-dissipation relation. In subsection [VC, we study the 
creation of particles by stochastic metric perturbations and see that this phenomenon can be related to the vacuum 
noise kernel. We show that the mean value of created particles is enhanced by the presence of metric fluctuations 
with respect to the same quantity in the "perturbed" semiclassical spacetime {M,gab+{hab)c)- 

Let us assume that the semiclassical spacetime (A^, ^ab) is stationary, i.e., that it possesses a global timelike Kilhng 
vector field i?^(/ab = 0, where £^ is the Lie derivative with respect to C"- Writing the Kilhng vector as C = {d/dt)"", 
this spacetime can be foliated by a family of Cauchy hypersurfaces St, labeled by the Killing time t, so we can give 
coordinates (t,x) to each spacetime point, where x are the space coordinates on each of these hypersurfaces. Using 
this foliation, we can construct a Hamiltonian operator H[g] in the way described in Appendix This is a time 
independent, i.e., independent of the Cauchy hypersurface Ej, Hamiltonian operator, so it represents the Hamiltonian 
operator in both the Heisenberg and the Schrodinger pictures. In this case, there is a natural Fock representation based 
on a decomposition of the field operator [g] in a complete set of modes of positive frequencies ujk with respect to C° , 
and their complex conjugates.^ This defines a natural Fock space, the many-particle states of which are eigenstates 
of the Hamiltonian H[g]. Thus, the notion of particles is physically well defined in this spacetime p],|60|,p[ . The 
Hamiltonian operator in this Fock representation, renormalized by normal ordering, is given by H[g] =X]fc'^fc (''k^k, 
where aj. and are the creation and annihilation operators on the Fock space. Here, the summation must be 
understood as representing either a sum over a set of discrete indices or an integral with some suitable measure (or 
a combination of these two possibilities). The time-evolution operator corresponding to this Hamiltonian operator is 
then given by U[g]{t,t')= e:iq){-iH[g] (t-t')). 

In this section, even if we sometimes write U or tf, we shall always consider these initial and final times in the limit 
ti— >— oo and tf—^+oo (we assume that such limits can be taken). 

A. The fluctuation-dissipation relation in a stationary background 

For a real scalar field quantized on the stationary spacetime (Al, 5ah), we can define a state of thermal equilibrium 
at temperature T. This state is described in the Heisenberg picture by the density operator of the grand canonical 
ensemble: 



^In some cases, additional restrictions may be necessary to avoid infrared divergencies, such as that the scalar field is massive, 
m^T^O, or that the norm of the Killing vector is not arbitrarily small |]l|,|6o|. 
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where /3 = l/ksT and ks is Boltzmann's constant (there are no chemical potential terms because we deal with 
a real scalar field). This kind of thermal states for fields in stationary curved backgrounds was first considered in 



Refs. 1^,^. Since the density operator (4.1) commutes with the time-evolution operator Z^[g](i, t'), the corresponding 
initial density operator in the Schrodingcr picture is simply p^{ti) = p[g]. 

Given any pair of operators in the Heisenberg picture, P[g]{x) and Q[g]{x), the expectation value {P{x)Q{x')) [g] 
depends on t and t' only through the difference t — t' , since 



{P{x)Q{x% 



Tr 



pP{x)Q{x'] 



Tr 



/5F^(x)e 



(4.2) 



where P^(x) and (3^(x) are the operators in the Schrodinger picture co rres ponding to P{x) and Q{x), respectively, and 
we use ( )^ to denote an expectation value in the state described by (4.1). In particular, with the choice p^{ti) — p[g], 

the kernels N^^'^''-[g]{x,x'), H^^'^'^[g]{x,x') and Hj^^'^[g]{x^x') depend on the time coordinates as a function of i — t' . 
Therefore, we can introduce Fourier transforms in the time coordinate as 



K{x,x') 



^e— 7^(c.;x,x'), 



(4.3) 



where K{x, x') is any function which depends on time only through t — t'. 

As it is shown in Appendix Wick's theorem can be generalized for th erm al A^-point functions, defined as 
expectation values of products of the field operator in the state described by (4.1). It is then easy to see that the 
expressions found in subsection piq also hold for the kernels Nf'''^\g]ix,x'), H^'^'^[g]{x, x') and Ef'^^jix^x') at 
finite T if we replace the Wightman and Feynman functions (3.38) by the analogous thermal expectation values. 

In this case, a simple relationship (in the for m of a fluctuation-dissipation relation) exists between the kernels 
N^'"^'^[g]{x,x') and H^^'''^[g]{x,x'). In fact, from (3.37), we can write these kernels as 



Nf^\x, x') = F^^^^x, x') + F^''^\x\ x), 8i H^'^^x, x') = ^''(x, x') - F^J^\x\ x), 



where we omit the functional dependence on gab- In terms of the Fourier transforms ( [4.3[ ), these relations are 

8 N:^^\oj- X, x') = F^^^^Lo- X, x') + F^J^\~Lo- x', x), 
8z H:1'\oj- X, x') = F^!^^^{lo; X, x') + F^„^'^\~io- x', x). 



(4.4) 



(4.5) 



By analytically continuing t to complex values in F^^'''^{x,x'), one can derive a symmetry relation for this bi-tensor 
which involves different values of this complex time. Taking into account that the time evolution of the operator tf:^ 
is given in this stationary case by i'^{t + At, x) = giHAt t<^{i^ x) e"*^'^*, and using the cyclic property of the trace, 
we get F^'"^'^{t, x; t', x') = F^'^°'^{t' , x'; i -I- i/3, x), or, equivalently, in terms of its Fourier transform. 



F^^iu;;^,^ 

This relation is known as the Kubo-Martin-Schwinger relation 
the following simple relation between N^^'^'^ and H^^^'^; 



(3lo 



^(-a;;x',x). 



(4.6) 



From this last expression and (4.5), we obtain 



(a;;x,x') = -i tanh ^ A^f ^'^(c.; x, x' 



which can also be written as 



/oo 
-OO 



abed {j.f! 



(t",x;t',x'), 



with 



duj 



sm{ujt) tanh 



^fcBTP[cosech(7rfcBTi)] , 



(4.7) 



(4.8) 



(4.9) 
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where P denotes a Cauchy principal value distribution. 

Since, as we have pointed out above, the kernels H^^'^ and N^'"''^ are free of ultraviolet divergencies in the limit 
n^4, we can define 



Hf'"^(x,x') = lim ^-2(n-4)^a6c<i/ ^^/N N'''"^'^ (x , x') = lim M"2(«-4)^afccrf. ^^/N 



(4.10) 



which are the kernels that appear in th e phy sical semiclassical Einstein-Langevin equation, Eq. ( P.llD, after performing 
the renormalization procedure in Eq. (^.3(: ). These physical kernels will also satisfy the relation ( tl.Sj) or, equivalently, 
their Fourier tra nsfo rms will satisfy (4.7). These results are independent of the regularization method used. 

The relation (4/7) can be written in an alternative way. Introducing a new kernel (this is actually a family of 
kernels) defined by H^^'"^{uj;k,x') = -iw 7^ '='*(w; x, x'), that is, H^^^'"^{x,x') = d-f^'""^{x,x')/dt, Eq. (|J) yields 



or, equivalently. 



Nf'"^{uj;x,x') = wcotanh 



I I 7n ('^,x,x ), 



(4.11) 



N^l'^^it, x; t', x') = / dt" J,^{t - t") ^if"'{t'\ x; t\ x'). 



(4.12) 



where 



Jo 



duj , f (3tjj 

— cos(cji) cotanh -— 
TT \ 2 



(4.13) 



This integral gives a distribution which is singular at t — Q and for t ^ reduces to J-p-D{t) = —-K [ksT cosech{TTkBT t)]^ 



The relations (4.7) or (4 



or the equivalent forms (4.11) or (4.12)] have the same form as the fluctuation-dissipation 
relations which appear in quite general models of quantum mechanics [p8|-^. The derivation of these relations is 
usually done in the framework of linear response theory, in which one considers the response of a quantum system, 
which is initially at thermal equilibrium, when an external classical time-dependent linear perturbation is "switched 
on." When evaluating the change in the expectation value of the relevant operator (the operator which couples to 
the perturbation) induced by the presence of the perturbation, a dissipative term can be identified as the term which 
changes the sign under a time reversal transformation in the perturbation. This term is characterized by a kernel 
called the dissipation kernel. It can be shown that the dissipation kernel is related to the fluctuations in equilibrium 
(in the absence of the perturbation) of the relevant operator by a relation which is exactly the same as (O) or 
(4.7). This is the fluctuation-dissipation relation. Using this linear response theory approach, the same fluctuation- 
dissipation relation has also been derived for some models of quantum many-body systems 43 1 or quantum fields 
p2| , ^ coupled to external classical fields. 

This fluctuation-dissipation relation appears also in the context of quantum Brownian motion (or "semiclassical" 
Brownian motion), in which one is interested in the dynamics of a macroscopic particle in interaction with a heat 
bath environment, usually modelized by an infinite set of quantum harmonic oscillators. In these models, when the 
variable representing the center of mass position of the macroscopic particle decoheres, it can be effectively described 
as a classical stochastic variable. The equation of motion for this stochastic variable is a linear Langevin equation 
with a Gaussian stochastic source. The classical variable introduced in linear response theory can be envisaged as the 
position of the Brownian particle, but now this variable becomes a dynamical stochastic variable. The dissipative term 
in this Langevin equation is the responsible for the irreversible dynamics of the Brownian particle. This term contains 
a dissipation kernel which is related to the correlator of the stochastic source by the relations (4.8) or (4.7) p^ , p^ . 
This is again the fluctuation-dissipation relation. There are also some models in which a purely quantum description 
of the Brownian particle is considered 1 63 6^ . The dynamics of this particle is then described by a quantum operator 
in the Heisenberg picture. By elimination of all the environment degrees of freedom in the equation of motion for this 
operator, one finds a quantum Langevin equation with quantu m f luct uatin g and dissipative terms. These terms are 
again related by a fluctuation-dissipation relation of the form (4.8) or (4.7). 



These analogies allo w us to identify the equivalent relations (4.8) and (4.7), and the analogous relations for the 
physical kernels ( 4.10| ), as the fluctuation-dissipation relation in our context. Because of this relation, the kernel 
H^'^'^{x,x') shall be called the dissipation kernel. The same fluctuation-dissipation relation was derived by Mottola 
J47[ in the context of quantum field theory in curved spacetime using the linear response theory approach. This 
author considered the case in which the background spacetime is static, but his result is easily generalized to a 
stationary background. In this paper, we have derived the same relation in the context of a Langevin equation for 
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stochastic metric perturbations, which would presumably describe the effective dynamics of gravitational fluctuations 
after a process of decoherence. For the particular case of a massless scalar field in a Minkowski background, this 
fluctuation-dissipation relation was derived in Refs. ||2^,Q from an explicit evaluation of the kernels. 

It is clear that the kernel N"-^'^'^{x,x') describes fluctuations in exactly the same sense as the quantum- mechanical 
models described above. In fact, as it was pointed out by Mottola from the point of view of linear response theory, 
it gives the fluctuations in equilibrium of the stress-energy operator. Alternatively, as we have shown in the previous 
sections, it gives the two-point correlation function of the Gaussian stochastic source in the semiclassical Einstein- 
Langevin equation. However, the term containing the "dissipation" kernel H^'"^'^{x,x') in the Einstein-Langevin 
equation does not generally change sign under a time-reversal transformation in the metric perturbations. 



1. Zero temperature limit 



A state of the scalar field which is of special interest is that described by p^{ti) — p[g] — |0)(0|, where |0) is the 
vacuum state. This vacuum state can be obtained as the zero temperature limit, T— >0, of the previ ous ther mal state. 
The fluctuation-dissipation relation for this state is easily obtained by setting T = in expression (4.7) or ( |4.8| ). We 
find x') = — «signw7V°^'^'*(u;;x, x'), or, equivalently, it has the form with 



duj 
2^ 



' sign uj = P 



(4.14) 



This fluctuation -dissi pation relation in the alternative form (4.11) reads A'^^'^'^(a;; x, x') = w signw 7°''^''(w; x, x'), or, 
it has the form ( 4.12 ), with 



(4.15) 



where Vf{l/t^) denotes a Hadamard finite part distribution, which is related to P(l/t) by 'Pf{l/t^) = —{d/dt)V{l/t) 
(the definitions of these distributions can be found in Refs. Bq]). 



2. High temperature limit 



Let us now consider the high temperature limit. This limit can only be performed when there exists a cutoff 
frequency fi, such that 7V^^'^''(u;; x, x') vanishes for ljj>VI (by ( |4.7| ), i/^^^'^'^(cj; x, x') will also vanish for these values 
of Lo). Such a cutoff frequency is usually related to a characteristic cutoff frequency of the environment degrees of 
freedom. The high temperature limit corresponds to the limit in which fc^T^Tifi. In this limit (keeping only the 
leading order contributions), we expect that thermal fluctuations dominate over quantum fluctuations. To study 
this limit, it is convenient to restore the dependence in Ti in the previous results. For this, o ne ha s to multiply the 
constants as and /3b by fi and the kernel H"^"^ by l/Ti in the Einstein-Langevin equation (3.36), and change the 
combination Plo by Tij3uj in the previous expressions. In this limit, we can approximate tanh(?i/3cij/2) ~ T%j3uj/2, and 
the fluctuation-dissipation relation reduces to 



^ Hf'"^iLo; X, x') = -i . 

;i ^ ' ' ' 2kBT 



■N. 



abed 



(w;x,x'). 



or, equivalently, since in this case Kpo{t) ~ {h/2kBT) (d/dt) S{t) 

1 d 



1 



i7f*"''(t,x;i',x') 



2kBT dt 



7V,f^''(i,x;t',x'). 



(4.16) 



(4.17) 



Note that {l/K)Hf'"^ is the kernel that appears in the Einstein-Langevin equation ( |3.36| ) when one writes the 
dependence in Ti explicitly. This relation has the same form as the classical Green-Kubo fluctuation-dissipation relation 
which appear s eith er in a classical theory of linear response 39 1 or in a classical theory of Brownian motion p^ , |67| . 
Notice, from (4.17), that in this high temperature hmit we can simply take ^^'^'^{x,x') — {h/2kBT) N^'"^'^{x, r). 
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B. The fluctuation-dissipation relation for conformal fields in a conformally stationary background 



In the case of a massless conformally coupled scalar field (m = and ^ — 1/6) and a conformally stationary 
solution of semiclassical gravity [for instance, a Robertson- Walker (RW) spacetime] , the fluctuation-dissipation relation 
derived in the previous subsection can be generalized when the state of the field in the background solution is the 



conformal vacuum or a thermal state built on the conformal vacuum. In this case, the action (3.1) for the scalar 
field is conformally invariant. It is convenient to preserve this conformal invariance when working in dimensional 
regularization. This can be done by changing in all the previous expressions which involve dimensional regularization 
the parameter ^ by the function g (n) = [n — 2)/[4(n — 1)] and, of course, taking m = 0. In this way, the dimensional 



regularized stress-tensor operator (3.6) is traceless. Let (A^,^^^) be a n dimensional conformally stationary spacetime, 
that is, a spacetime with a global timelike conformal Killing vector field C^: £jjab — (2/?^) V'^Cc 5a6j where Va is the 
covariant derivative associated to g^^. This means that the metric g^f, is conformally related to a stationary metric 
9ah'- Vabi^) — ^^^'"^''9ab{x), whcrc zu{x) is a scalar function. As previously, writing — (d/dt)"', the semiclassical 
spacetime can be foliated by Cauchy hypersurfaces St and coordinates {t, x) can be assigned to the spacetime points. 

There is a "natural" Fock representation based on a decomposition of the field operator ^nig] in terms of a 
complete set of modes {ukj^x)}, solution of the Klein-Gordon equation with metric 'g^^,, of the form Uk^J^x) — 
g-(ra-2)t37(x)/2y^^^^^^ ^j^gj,g jy^^^j;)} is a complete set of mode solutions of the Klein-Gordon equation in (M.gab) 
which are of positive frequencies uj^ with respect to Hence, in this sense, we can write the field operator as 
$„[g] — e~'^"~^^™/^4>„[(7], where ^nlg] is the field operator in the stationary spacetime {Ai,gab)- Assuming that no 
infrared divergencies are present, so that this quantum field theory construction is well defined, the conformal vacuum 
|0) is defined as the vacuum state of the Fock space corresponding to this representation. If and fifc are the creation 
and annihilation operators on this Fock space, this state satisfies dk\0) = 0. As shown in Appendix in this case 
we can construct a conserved energy operator which can be identified with the Hamiltonian of a field quantized on 
{M,gab)- E[g] ^ H[g] ^J2k^k o.l.ak. This energy operator, however, is not a time-evolution generator for the field 
operator $„ [g] , it generates the time-evolution of the conformally related operator $„ [g] . The many-particle states 
of the Fock space built on the conformal vacuum are eigenstates of this energy operator. 

From this energy operator, a state of thermal equilibrium for a conformal scalar field quantized on (A^,5af,) can 



be defined using the density operator (4.1). Thermal equilibrium states defined in this way were first proposed by 
Gibbons and Perry |Q . These authors were inspired in a result by Israel |Q in the framework of relativistic kinetic 
theory, who found that thermal equilibrium distribution functions can be defined for massless particles in conformally 
stationary spacetimes. A number of applications have been developed in the literature to study finite-temperature 
effects of quantum conformal fields in RW universes ||6^] and in two-dimensional spacetimes [ [70[ . 

Let us begin with a solution of the semiclassical Einstein equation ( |3.10| ) con sisting of a quantum conformal scalar 
field in a conformally stationary spacetime (Al,g^j), in the thermal state ( [4.l[) . Taking into account that the action 



(3.3) with m — and ^ — ^(n) satisfies Sm[g,'^n] = Srn[g,^n], where = e ^)^/^<i>„, it is easy to see that 
'^n''[g] = e~'"+^'™r°^[(7]. Therefore, the kernels evaluated in the thermal state at a temperature T can be related to 
the corresponding kernels for the stationary background {A4,gab)- For the noise kernel, we have 

7Vf^'^[g](x,x') = e-("+2)^(")e-("+2)^("') A^f^'^[g](x,x'), (4.18) 
and the same relation holds for the kernels H^^'^'^ and H^^"^- Since the kernels iV^'"^''[5] and H^^'^[g] satisfy the 



relation ( [4.8| ) [or, equivalently, (4.7)], this leads to a fluctuation-dissipation relation between the kernels N '^ \g \ and 



H^^'^[g]. The same relation holds for the physical kernels obtained by taking the limit n-^ 4 as in Eq. ( |4.1C ). For 



the conformal vacuum state, which corresponds to T = 0, the fluctuation-dissipation relation follows directly from 



the result of subsection IV A 1 . In the particular case of a spatially flat RW solution of semiclassical gravity, this 
conformal vacuum fluctuation-dissipation relation was obtained before in Ref. p9| ] after an explicit calculation of the 
corresponding kernels. The same relation was derived in Ref. [ [T^ in the framework of a "reduced" version of the 
Einstein-Langevin equation inspired in a Bianchi-I type "mini-superspace" model. 



C. Particle creation 



Let us now return to the case in which (A^,(7ab) is stationary, the scalar field has arbitrary mass m and arbitrary 
coupling parameter ^, and consider the stochastic perturbation hab- Note that (Al,5ab + ^ab) can be viewed as 
representing an ensemble of spacetimes distributed according to some probability distribution functional. We are in 
fact considering a scalar field quantized on each of these spacetimes, described by the operator ^[g+h], and the family 
of states of the field, described by p[g + h]. 
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Let /iab[C] be a solution to the semiclassical Einstein-Langevin equation, Eq. (2.11), whose moments vanish for times 
t < ti or, at least, they vanish "asymptotically" in the remote past (t— >— oo). This means that there is a "remote 
past epoch" {t < ti or t ^ —oo) in which hab behaves deterministically as a zero tensor. In that case, if we take 
p^{ti) = |0)(0|, where |0) is the vacuum of the natural Fock space for the field quantized on {Ai,gab), and we consider 
the limit ti —^ — oo, we have p[g + h] — |0, in) (0, in|, where |0,in) represents the family of "in" vacua for the field 
quantized on {M, gat + hat)- Treating hat as a classical "external" perturbation, one could construct a Hamiltonian 
operator H[g + h]{t) in the Heisenberg picture for which |0,in) would be the ground state in the "remote past epoch." 
However, at later times, due to the presence of the perturbation hat, this "in" vacuum state will generally not be the 
ground state of the Hamiltonian. One then says that "particles" are created in the "in" vacuum. 

Physically meaningful many particle "out" states, in particular, an "out" vacuum |0, out) for the scalar field in each 
of the spacetimes {M, gab+hab) can be defined if there is also a "far future epoch" for which hab vanishes (in the same 
statistical sense as above), either in an exact way for t > tp oi "asymptotically" for t^+oo. When this is the case, 
the vacuum persistence amplitude (0, out|0, in) [g + ft,] = e'^[^+'*l is given by the following path integral: 



y"l?[$„] (0,out|$„(t2),t2)[g+/i] ($„(ti),ti|0,in)[.g + /i] e^^^ [<,+/.,*„] ^ (4 ^g) 



where \tpn, ii) and |(p„, denote, respectively, eigenstates of the field operator $„[g+/i](i, x) at some arbitrary times 
t = ti and t — t2, where t2 > ti, with eigenvalue (p„(x), and where the integration domain for the action is between 
ti and t2- The wave f unctio nals (iy9„, ii|0, in) and (</?„, i2|0, out) have in general a dependence on the metric, which 



we have indicated in (4.19). In the limit ti —00 and t2 +00, these wave functionals do not depend on the 



perturbation hab- The total probability of particle creation is given by |71| 

P[h;g] = 2 lim ImW[g + h]. (4.20) 



+4 



One can show that ImW^ is free of ultraviolet divergencies in the limit n— >4, and that it is always positive or zero, 
so that the probability P is well defined by this expression. 

As we have done in the previous section for the influence action, we can now expand the action ^^[^ + /i] in the 
perturbation hab- In order to do so, one has to evaluate the functional derivatives of + h] in the background 



metric gab- Using (4.19), these derivatives can be related to "in-out" matrix elements of operators in the background. 
Since gab is stationary, the "in" and "out" vacua in the background must be identified with the natural vacuum |0). 
Therefore, these background "in-out" matrix elements become expectation values in the state |0). It is then easy to 
see that the expansion of M^[(7-t- /i] in the metric perturbation hab is equal to that of Sip[g + h'^ , g + h~] with /i^^ = hab 
and h~^ = 0, and taking the expectation values in |0). In particular, from the imaginary part of this expansion [see 



Eq. (3.18)], we get 



P[h-g]= d^xdSV^9(x)V^^)hab{x)N'^'''''[g]{x,y)h,d{v)+0{h''), (4.21) 



where ]\[°-'"^'^ is the zero temperature physical noise kernel define d in ( 4.10 ). This physical noise kernel is related to the 
lowest order quantum stress-energy fiuctuations in vacuum by ( |2.9| ). Not e that the higher order corrections in (4.21) 
would contain higher order vacuum stress-energy fluctuations. Eq. ( [4.21 ) is a generalization of an expression derived 



by Sexl and Urbantke | |72| for the total probability of particle creation by metric perturbations around Minkowski 
spacetimc. 



Eq. (4.21) gives also the expectation value of the number operator for "out" particles in the "in" vacuum, computed 
to lowest order in the metric perturbation. In order to show this, let us expand the scalar field action as the action 
in the stationary background plus interaction terms (the terms containing the metric perturbation) . The interaction 

term to lowest order in hab is S^^'' — J d^x Cn'^ [$„, h; 5], with Cn^ — (1/2) ^—gT°-^[g, <&„] hab- In order to construct the 
5-matrix operator, we need the interaction Hamiltonian density operator in the interaction picture. Note that the field 
and canonical momentum operators in the interaction picture can be identified with the operators $n [9] and H„[5], 
respectively. Following Appendix we can obtain the canonical Hamiltonian density for the metric gab+Kib and work 
out the interaction term to first order in the metric perturbation. Although in this case the interaction Lagrangian 
density depends on the derivatives of the field, we find that, to first order in hab, the interaction Hamiltonian density 
operator in the interaction picture is given by — £l^''[$„[g], h; g]. Hence, to first order in the metric perturbation, the 
S'-matrix operator is given by 5 = 1 + S^^^ +0{h'^), where 

i 



S^^^ = - I d^xV^f,f[g]hab- (4.22) 



18 



The expectation value of the "out" particle number operator, 7V°"', in the "in" vacuum (in the Heisenberg picture) 
is given by N[h;g] = (0, in| iV°"* |0, in) = (0|S'WS'|0), where N is the particle number operator in the background 
N = J2k ^fcOfe- To lowest order, we have 



iV[/i;g] =^|(U,lp|5(i)|0) +0{h'), 



where Ip) is the two-particle state |lfc, 1^ 



k,p 



(4.23) 



|0). Clearly, since 5^^' is quadratic in the field operator, at this 



order N can also be written as N/2 = T^^jOlS^^^^n) {n\S'-^^ \0) - {0\S^^^^\0){0\S^^^\0) +0{h^), where {\n)} represents 
the complete orthonormal basis of the F ock s pace. Using ( |4.22| ), this last expression can be written in terms of the 
vacuum noise kernel Nf'"^[g]{x,y) [see (| .14|) ] . Taking the limit n — > 4, we see that the expression for one half of 
the number of created particles N[h;g]/2 to lowest order in the metric perturbation coincides with that for P[h; g] in 
Eq. ( ^^ . 

The energy of the created particles, defined as E[h; g] = (0, in| J2k ^kN\ 
"out" number operator in the k mode and H = ^y. ivk a^ak is the Hamiltonian operator in the background, is similarly 
given by 



^out 



0,in) = (0|5tiJ5|0), where iV°^* is the 



k,p 



(4.24) 



Comparison of ( [4.24| ) with (4.23) and ( [f.2l| ), suggests that it may be possible in some cases to write this last expression 
(in the limit n— >4) in terms of the Fourier transform of the vacuum noise kernel. 

As an example, let us consider the case when {M.,gab) is (M"*, r/afc) [^llj73| , which is the trivial solution of semiclassical 
gravity. Working in a global inertial coordinate system {x^}, in this case the kernels depend only on the difference 



ip-(x-y) 



K{p), where p-x \ 



{x—y)^ and, thus, we can define their Fourier transforms as K{x—y) = (27r) d'^pe 
Introducing the Fourier transform of hab{x) in a similar way [note that /iaf,(— p) = /i*j(p)], ( 4.21 ) can be written as 



N'^'^\p)hUp)Ka{p) + 0{h' 



On the other hand, the energy of the created particles is given by |74| 

d^p 



{2^Y 



(4.25) 



(4.26) 



The vacuum noise and dissipation kernels for a Minkowski background can be written in terms of two pairs of scalar 
kernels, Nr{x — y) and Dr{x — y), respectively, with r — 1,2 [|7^ (see also Ref. [^ for a particular case i n which 
7V2 = D2 = 0). Each pair of kernels {Nr, Dr) satisfies the fluctuation-dissipation relation found in subsection IV A 1. 
One finds [|l]j74) that 



R''>ip) 'n2{p), 



(4.27) 



where Cabcdip)' R^'''{p) and Nr{p) are, respectively, the Fourier transforms of the linearized Weyl tensor, the scalar 



curvature and the kernels Nr{x — y), r - 
fluctuation-dissipation relation, that 



1,2; Nr{p) depend only on p 



■'hvP'^p" 



It is then easy to see, using the 



E[h-rf, 



dp 
(27r)4 ^ 



R''\p) ^Mp) 



(4.28) 



Hence, in the case of a Minkowski background, the energy of the created particles can be expressed in terms of 
the dissipation kernels Di and D2 for the Minkowskian vacuum. It is not clear, however, that, for other stationary 
backgrounds, the energy of the create d par ticles can be related to dissipation in vacuum in a similar way. 

The probability of particle creation (4.21) is a fluctuating quantity, due to the functional dependence on the stochas- 
tic p ertur bation hab- We may compute its averaged value {P[h]g])c, which [neglecting the higher order corrections in 
Eq. (4.21)] is given by 



P[h;g]), = P[{h),-g] + / d'x d'y y/^^^/^) N'^'^^%]ix,y) {hUx)hUv))c 



(4.29) 
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where /i^^ = hab — {hab)c- The first term in the right hand side of Eq. ( 4.29| ) is the probabihty of particle creation (or 



one half of the number of created particles) that one would obtain in the spacetime {M, gab+{hab)c)- The second term 
will be greater than zero when stress-energy fluctuations are present^ since, from the Einstein-Langevin equation, this 
implies (/iaft(a;)ft.^^(y))c 7^ 0. Note that, when this is the case, from the fluctuation-dissipation relation of subsection 



EVA 1, the vacuum dissipation kernel will be also non- vanishing. Hence, metric fluctuations induced by matter stress- 
energy fluctuations generally increase the mean value of the number of created particles with respect to the same 
quantity in the "perturbed" semiclassical spacetime (Al, (7ab+ (/ia6)c)- 

The above result for the total probability of particle creation and number of created particles can be easily general- 
ized to the case of a massless conformally coupled scalar field and a conformally stationary semiclassical background. 
When this background is a spatially flat RW universe , performing conformal transformations in the metric pertur- 
bations and in the kernels as in Eq. (4.18), one gets expressions analogous to (4.25), (4.27) and ( 4.2^ ) with N2 = D2=0 



(see Refs. [p6Hl9|j3l| for more details) 



V. CONCLUSIONS 



In the first part of this paper, we have shown how a consistent stochastic semiclassical theory of gravity can be 
formulated. This theory is a perturbative generalization of semiclassical gravity which describes the back reaction of 
the lowest order stress-energy fluctuations of quantum matter fields on the gravitational field through the semiclassical 
Einstein-Langevin equation. We have shown that this equation can be formally derived with a method based on the 
influence functional of Feynman and Vernon, where one considers the metric field as the "system" of interest and 
the matter fields as part of its "environment" ^ . Our approach clarifies the physical meaning of the semiclassical 
Langevin-type equations previously derived with the same functional method |p^- |l^ , |l9| -p5| , since it links the source of 
stochastic fluctuations to quantum matter stress-energy fluctuations, and allows to formulate the theory in a general 
way. At the same time, we have also developed a method to compute the semiclassical Einstein-Langevin equation 
using dimensional regularization. This provides an alternative and more direct way of computing the equation with 
resp ect to the previous calculations, based on a specific evaluation of the effective action of Feynman and Vernon 
|p^- |l4| , p^ -|2^] . In a subsequent paper Q , we shall apply this method to solve the Einstein-Langevin equation around 
some simple solutions of semiclassical gravity. 

The second part of the paper was devoted to the existence of fluctuation-dissipation relations and to particle creation 
in the context of stochastic semiclassical gravity. When the background solution of semiclassical gravity consists of 
a stationary spacetime and a scalar field in a thermal equilibrium state, we have identified a dissipation kernel in 
the Einstein-Langevin equation which is related to the noise kernel by a fluctuation-dissipation relation. The same 
relation was previously derived by Mottola |47| using a linear response theory approach. We have also generalized this 
result to the case of a conformal scalar field in a conformally stationary background solution of semiclassical gravity. 

Our analysis seems to indicate that for a fluctuation-dissipation relation to be present in stochastic semiclassical 
gravity, the semiclassical background solution must satisfy certain conditions. In this paper we have just analyzed the 
simplest cases for which such a relation exists. Further work must be done to investigate whether a similar relation 
is present in other situations of physical interest, such as black hole backgrounds ^8 49, 2^], or non-conformal fields 



in RW backgrounds in the instantaneous vacua or the thermal states defined in Ref. |75[|. 

We have also studied particle creation by stochastic metric perturbations in stationary and conformally stationary 
(for conformal matter fields in this latter case) background solutions of semiclassical gravity. We have expressed 
the total probability of particle creation and the number of created particles (the expectation value of the number 
operator for "out" particles in the "in" vacuum) in terms of the vacuum noise kernel. We have shown that the 
averaged value of those quantities is enhanced by the presence of stochastic metric fluctuations. In the particular 
cases of a Minkowski background and a conformal field in a spatially flat RW background, the energy of the created 
particles can be expressed in terms of the vacuum dissipation kernels. 

It should be stressed that the concept of particle creation is only well defined when the solutions of the Einstein- 
Langevin equation vanish in the "remote past" and in the "far future" (at least, "asymptotically"). However, there 
can be physically meaningful solutions of the Einstein-Langevin equation that do not satisfy these rather strong 
conditions. In this case, vacuum noise and dissipation in stochastic semiclassical gravity can include effects that are 
not associated to particle creation. 



^Except in some rare cases, for which N'^'"^''{x,y) is not strictly positive definite and {h^ab{^)^\d{y)) ^ is such that it "hits" the 
zero eigenvalue. 
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APPENDIX A: EXPANSION OF THE STRESS-ENERGY TENSOR AROUND A BACKGROUND 

METRIC 

The expansion of the stress-energy tensor functional 

around a background metric gab is given by 

T'^'ig + h, $„] = T-'[g, $„] + r'^^'ig, h] + 0{h'), with 

T(^'"^[g, $„; h] = -T-'[g, $„] - T^^lg, $„] (v'I'n Vc<f„ + m'^l) h^" + 1 g'^' V^'&n v'*n 

+ (v"/ic + V'K - Vch'"' - 23^" v'^crf + g"' Vch',) +2/^"" □ - 2g^'^Kd V 

where the covariant derivatives and curvature tensors are those of the metric gab, and indices are raised with inverse 
background metric 5"''. 

APPENDIX B: HAMILTONIAN OPERATOR IN A STATIONARY SPACETIME AND THERMAL 

WICK'S THEOREM 

In this appendix, we construct the Hamiltonian or energy operator for a quantum scalar field in a stationary 
spacetime. For a more rigorous mathematical treatment, see Ref. [pOj and, for the particular case of a static spacetime, 
see Ref. We also show how this construction can be generalized for a conformal scalar field in a conformally 
stationary spacetime. Using this Hamiltonian to define a thermal density operator, we shall see how thermal four- 
point functions can be expressed in terms of thermal two-point functions ("thermal Wick's theorem"). 

Let {M,gab) be a n dimensional stationary spacetime, that is, a spacetime with a global timelike Killing vector 
field = (d/dt)'^, and consider a linear real scalar field <!>„ on it. Assuming that the spacetime is globally hyperbolic, 
we can foliate it by a family of Cauchy hypersurfaces Sj, labeled by the Killing time t (hypersurfaces of constant 
t), and give coordinates to each point of the spacetime x'^ = (t, x), where x = (a;*) are local coordinates on each of 
these hypersurfaces. Let be the future directed unit {i.e., Uan"" = — 1 and n* > 0) vector field normal to each 
hypersurface Ej. The induced metric on each Et by the spacetime metric is Qab = gab + naUi, |^^, then is a 
projector orthogonal to n". We can decompose the Killing vector into its normal and tangential parts to each Ej: 
C = Nn'^ + N°-, where N = —C^ria and N'^ = are, respectively, the lapse function and the shift vector. In 

the basis associated to the coordinate system {a;''}, the components g^^i, of the metric are independent of t and can 
be written as gu — —N"^ + NiN\ gu = Ni, gij = qij, with Ni = qijN^ . One can also write \/—g = where 
q = det (qij). 

To construct the classical Hamiltonian, we write the Lagrangian density as 

= I [{nf^dM^ - q''d,<i>ndj<i>n - (m^ + ^R)^l] , (Bl) 

where n* = = —N^/N and g*^ is the inverse of qij, qikq^-' — S^. The momentum canonical conjugate to 

$„ is H„ = y/q n^dfi'^n and the Hamiltonian density is constructed as usual, 7i„ — H„ — £„, from which the 
Hamiltonian functional on the hypersurface Et is given by Hit) — d^^^x-Hnix). Integrating by parts and dropping 
surface terms, we get 
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7? 




(B2) 



where Di is the covariant derivative on the n—1 dimensional Riemannian spaces {T,t,qij) (associated to the metric 
Qij), and A = D'^Di is the associated Laplace-Beltrami operator. For a field and its momentum conjugate n„ 
satisfying the Hamiltonian equations of motion, this Hamiltonian is a conserved quantity, i.e., it is independent of t. 

The analogous quantity in the Lagrangian formalism is the canonical energy functional, which is defined in terms 
of the canonical stress-energy tensor functional. 



Tarla, = Va*n Vb*« - J 9ab V^*n Vc*n ^ J 9ab (jn^ + ^R) 



(B3) 



as Ec. 



dT, vf^C,^ T^b^[g, $„], where E is a Cauchy hypersurface, rf" is the future directed unit vector field normal 



to S, and dS is the invariant volume element on S constructed with the metric induced by gab. By Noether's theorem 
this energy functional is conserved, i.e., it is independent of the choice of Cauchy hypersurface S, when $„ 
satisfies the Klein-Gordon equation. Choosing E( as the Cauchy hypersurface, we can obtain an expression for Ecan 
after the substitution of n„ by ,/q n'^d^^n h^ the Hamiltonian (B2). Note that we can also introduce an energy 
functional E = J.^dT, n°'('' Tab[g,^n], where Tab is the stress-energy tensor functional (2.2) ||5^. For a field $„ 



satisfying the Klein-Gordon equation, this is also a conserved quantity. However, choosing a Cauchy hypersurface Et, 
one can show that n°'('' {Tab ~ 2^afc") ^ divergence on the space (E^, g^) and, thus, dropping surface terms, we have 

E — Ecan. 

We can now formally construct the Hamiltonian "operator" in the Heisenberg picture simply by replacing $„ and 
H„ by their corresponding operators $„ and n„ in (B2) and using, as always, a Weyl ordering prescription for the 
operators. This operator is a conserved quantity, that is, it is independent of the time t; therefore, it is equal to the 
Hamiltonian operator in the Schrodinger picture and we simply denote it by H. Since the momentum operator in 
the Heisenberg picture satisfies H„ = this Hamiltonian operator can also be obtained from the canonical 

energy functional [hence, H represents also a conserved energy operator]. Taking into account that the field operator 
satisfies the Klein-Gordon equation, we find 



H = - f d"-^xVg^ 



St 



(B4) 



In this case, there exists a natural Fock representation based on a decomposition of the field operator $„ in terms of 
a complete set of modes {ufe^(a;)}, solution of the Klein-Gordon equation, which have positive frequency with respect 
to the Killing vector = (9/9t)°: dtUkJ^x) = —iujkUk^(x), with ujk>0. The label k of each mode must in general be 
understood as representing a set of discrete or continuous indices, and, thus, the summations over k represent either a 
discrete sum or an integral with some suitable measure (or a combination of these two possibilities). We assume that 
these modes have the same physical dimensions as the field (this is the reason why we put a subindex n). These 



modes have to be orthonormal with respect to the inner product 



2) = ~* /s ' 



hda' 



4 



'da 



is independent of the Cauchy hypersurface E when 0i and 02 are solutions of the Klein-Gordon equation 
(Mfe„ = ^ki and (ufc„ , uf^ ) = 0. 

The field operator can then be written as 



j, which 
0, 



I.e.. 



E 



Ukix)dk+ul(x)d^ 



(B5) 



where dt and dk are creation and annihilation operators on the Fock space associated to this mode decomposition. 



which satisfy the usual commutation relations [ffUTq]. Using th ese commutation relations and the orthonormality 
conditions for the modes evaluated on E(, substituting (B5) into (B4), one finds the Fock space representation of the 
formal Hamiltonian "operator" H = '^^u)k{d\dk -\- \). We can make this last expression well defined by subtraction 
of the "divergent" constant c- number X)fc(^fe/2)i that is, we can introduce a renormalized Hamiltonian operator as 



E 



(B6) 
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Note that Hji is given by an expression similar to (B4), but adding a normal ordering prescription for the operators 
Ofc and a\, or, equivalently (dropping surface terms), by Hr = J-^dY, n^C'' : Tn^^^ [g] ■ — f -""-^^ 
fn^f\g\ is defined in (3^), r^^^^iff] is analogously defined after (^3|), and 



n-C'':T-;>[5]:, where 



means normal ordering ||76|,|5| . The vacuum 
and the many-particle states of the Fock space are eigenstates of this Hamiltonian operator with zero and positive 
eigenvalues, respectively, (given by the sum of the w^'s corresponding to the particle contents of the state). 
From (B5) and Hn = ^/qn^d^^^m using the positive frequency condition and 



R,ak\ 



i^k Ofc, 



(B7) 



we get dt^n — i [Hr, $„] and dttln [Hu, tin]- These are the Hcisenberg equations of motion, which are equivalent 
to the Klein-Gordon equation for the operator $„. From these equations, we see that the operator exp {~iH]i{t — t')) 
generates the time evolution of operators in the Heisenberg picture. 



1. Conformal field in a conformally stationary spacetime 

Let us now consider a massless conformally coupled real scalar field $„ in a n dimensional spacetime {A4,'Qj^ 



assumed to be conformally stationary and globally hyperbolic. The action Sm[g,^n] for the field is given by (3.3) 
with m = and £^—^{n) = {n — 2)/[4(n— 1)]. In this case, the spacetime has a global timelike conformal Killing vector 
field C° = {d/dtY, which satisfies VaCb +V&Ca — {'^/''T')\/%c'gabi where Va is the covariant derivative associated to the 
metric g^^i,. The metric g„,h is conformally related to a stationary metric gab- Vabi^) = e^^'-^-'(7ah(a;). The foliation of 
the spacetime by Cauchy hypersurfaces Et and the coordinates a;'' = (i, x) are introduced as above. 

Given a Cauchy hypersurface E, with unit normal (as above, we take n* > 0), we can introduce the energy 
functional as £^ = J.^ dTj n'^''Q^ Tab[gT^n\, where dE is the invariant volume element constructed with the metric on 
E induced by the metric 'g^^. Given that the stress-energy tensor Tab[g, is traceless when the field satisfies 
the Klein-Gordon equation, it is easy to see from the equation for the conformal Killing vector C,"' that this energy 
functional is conserved. In fact, choosing a hypersurface Et to evaluate this energy, and introducing <i>„ = e("~^)'^/^ $„, 
it is easy to see Q that 

[ rf"-ixy^n'^C'T,b[g,<i>„], (B8) 

where n° and qij are constructed with the metric gab- Thus, E is equal to the energy functional for the field $„ in 
the stationary spacetime {M.,gab)- 

Using the "natural" Fock representation, based on the decomposition of the field operator in terms of modes 

UkJ^x) — e^'"^^^^('^)/^Ufc^(a;), we can construct the renormahzed energy operator in the Heisenberg picture, ER[g\^ 
associated to the above energy functional E. Here, as above, {ukjyx)} is a complete set of modes, solution of the Klein- 
Gordon equation in the stationary spacetime {AA,gab), which are of positive frequency with respect to = {d/dt)°'- 
Dropping surface terms, we get Enlg] = J2k (^iflk-, where a[. and Ofc are creation and annihilation operators on the 
Fock space associated to these conformal modes. 

Alternatively, one can perform the transformation $„ = e^*^"^^^"/^ $„ in the scalar field action, which is then 
transformed to 5'm[5, 'I'nlj and construct the Hamiltonian associated to this transformed action, which is given by the 
above expressions for the stationary case. This is equivalent to making a canonical transformation in the Hamiltonian 
formulation of the theory. One then introduces an operator ^n\g\ = e'^"^^)™/^$„[g], which can be identified as the 
field operator in the Heisenberg picture in the stationary spacetime {M.^gab)- The associated Hamiltonian operator 
can be identified with the operator iif^g\ constructed above (and, obviously, it coincides with Eb\§^- Note that this 
Hamiltonian or energy operator generates the time evolution of the operator \g\ rather than that of the "physical" 



field A generalization of this last approach has been used in Ref. |75| for scalar fields with arbitrary mass 

and arbitrary coupling to the curvature in a RW spacetime to construct a time-dependent Hamiltonian operator 
whose ground state at each fixed instant of time is a Hadamard state. A similar construction starting with the above 
energy functional is given in Ref. In the massless conformally coupled case, these time-dependent Hamiltonian 
constructions reduce to the construction sketched in this appendix. 
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2. Wick's theorem for thermal states 



From the Hamiltonian op erator ( ]B6| ) (here, we drop the subindex R), we can define a state of thermal equihbrium 
for the scalar field as in (4J). Following partially the proof presented in the appendix of Ref. [Q, we shall ne xt sh ow 
how Wick's theorem can be generalized for the associated thermal TV-point functions. First, note that, from (B7), 



where a — +,—, a"^^ =ak and a\. ^ = a\. Using this and the cyclic property of the trace, we get 



'(7) 



(B9) 



(BIO) 



where we have used the commutator [a'^\ O;^^] to represent either 5ki, —Ski or (such commutator does not represent 
an operator in the last equation) . Writing the field operator $„ (x) in terms of the operators a[,"^ , the associated four- 
point thermal functions can be expressed in terms of (a^"''a['^''ar*''as'^'')r • Taking into account that the commutator 
[a^k'\ o,['^'^] is a c-number, one has the following identity 



-(a).(7).(5).(o-) 



.{a) ,(7) 



On the other hand, from (B9) and the cyclic property of the trace, we have 



(Bll) 



(B12) 



Using the last two equations, we get 



1 



(1 _ g-a/3wfc^ 



-(7) 



(B13) 



which, from (BIO), yields 



+ (ai")a(^))Ja|-)a(^))^, (B14) 
and, hence, we have 

[^nixi)^n{x2)^n{x3)^n{x4:)) = (a^l )*n (2:2 ) ) ($„(a;3)'^n(a;4)) + (a^l )$n (a^s) ) ($n(a;2)<l„(2;4)^ 

+ (^„(a;i)^n(x4)\ (4„(a;2)4„(x3)\ . (B15) 



A similar expression holds for the four-point function of time-ordered products (consider ti > ^2 > ^3 > ^4 in the last 
equation). These results can be easily generalized to thermal 2A^-point functions (A^elsl). On the other hand, from 



(B9), we can see that {al^'')r^=0 and, following similar steps, we can show that the thermal (2iV— l)-point hmctions 
vanish. 
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